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Decide_param(T, R){ 

Nb=i; 

for(Np=2;Np<=Npmax;Np=Np*2) 

if(l/R <Tmin+TlQ+TlDCT 

+(Tc+TM-Pmax)(Np-l)/NpNb+TACDCmax){ 

Np=Np/2; 
return(Np,Nb); 

} 

Np=Npmax> 

for(Nb=2;Nb<=Nbinax;Nb=Nb*2){ 
Np=Np/2; 

ifCl/R <Tmin+TlQ+TlDCT 
+(Tc+TM-Pmax)(Np- l)/NpNb 
+(Tc+TM-Bmax)(Nb- l)/Nb+TACDCmax){ 

Nb=Nb/2; 
return(Np,Nb); 

} 

} 

Nb=Nbmax; 
Np=Npmax/Nbmaxj 

return(Np,Nb); 
} 
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